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Abstract. We examine the development of a transient 
spiral arm in a disk galaxy made up of both gas and 
stars. To this end we have performed numerical simula- 
tions in a shearing sheet (basically a rectangular patch of 
a disc) that contains gas in the form of clouds behaving 
like Brahic's (1977) sticky particles, and stars that ap- 
pear as a background continuum providing the perturba- 
tion forces. These are computed from the theory of swing 
amplification, using Fuchs' (1991) work. We describe the 
evolution of our model under a single and under recur- 
ring swing amplification events, discerning three phases. 
Furthermore, we give an interpretation of this evolution 
in terms of a variation of the epicyclic frequency with the 
distance to the wave crest. We also assess the importance 
of self gravity in the gas for our results. 

Key words: Galaxies: kinematics and dynamics ~ ISM: 
kinematics and dynamics 



1. Introduction 

Since the work of Orr (1907) it has been known that cer- 
tain perturbations in stable shearing flows show a tran- 
sient exponential growth. This phenomenon was discussed 
extensively by Goldreich & Lynden-Bell (1965) and Julian 
& Toomre (1966) in the context of galactic dynamics, us- 
ing hydrodynamics and statistical mechanics, respectively. 
Both papers examine the dynamics of a small patch within 
a larger disk and find strong transient growth of density 
perturbations as they "swing by" from leading to trailing. 
Toomre (1990) has argued that this amplification mech- 
anism is the principal dynamical process responsible for 
spiral arms in disk galaxies. 

While gas is known to play a major role in quasi- 
stationary density waves (Roberts 1969), there has been 
little research on the dynamical behavior of the interstel- 
lar medium (ISM) in transient density waves. It is this 
issue the present paper addresses. The key question we 
want to answer is whether the short lifetime of a swing- 
amplified perturbation still allows a noticeable response of 
the ISM. Additionally, we look for kincmatical signatures 



that might allow a discrimination of quasi-stationary ver- 
sus transient density waves. 

So far, the only study of swing amplification in a two- 
component medium was undertaken by Jog (1992) who 
extended the formalism of Goldreich & Tremaine (1978) to 
include two fluids with different stability numbers. While 
these works used Lagrangian coordinates to solve the dy- 
namical equations, in this paper we employ a different 
approach using Eulerian coordinates. This enables us to 
consider a succession of swing amplification events in a 
consistent way. 

The outline of this paper is as follows: In section 2 
we describe our model with a brief review of the results 
of Fuchs (1991) important to this work and a discussion 
of our treatment of an ISM consisting of discrete clouds. 
Section 3 presents the main results. In a first part, we 
describe the behavior of our model under a single pertur- 
bation. This behavior is interpreted in the second part. 
After a discussion of the model with repeated perturba- 
tions, section 3 closes with an investigation of the effects of 
self gravity in the ISM. In Section 4 we briefiy summarize 
the main conclusions we draw from this work. 



2. The Model 

Due to the inherently local nature of swing amplification 
no global disk model is required. Instead, we use a sliding 
box scheme similar to the one used by Toomre & Kalnajs 
(1991). The central idea of the sliding box scheme is to 
examine the dynamics of a rectangular patch of a galac- 
tic disk that is periodically repeated such that when, for 
example, a particle leaves the patch at the rotationally 
leading border, it is re- fed at the trailing one. This resem- 
bles the periodic boundary conditions used in plasma and 
solid state physics, with the complication that in galactic 
dynamics the periodic continuations must share the gen- 
eral shear of the disk. In consequence, particles leaving 
the patch radially must be re-fed at the opposite edge at 
positions varying with time, and the difference in circu- 
lar velocities between the inner and outer edges must be 
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accounted for, so as to keep the accelerations continuous 
during the re-feeding process. 

In such a patch pseudo-Cartesian coordinates x — 
r — Tq and y ~ ro(?? — ^oO ^'"^ introduced, where r, d 
are galactocentric polar coordinates, ro is the distance of 
the center of the patch from the center of the galaxy, flo is 
the angular velocity at rg, and t denotes the time. After in- 
serting these new coordinates into the equations of motion 
derived from the Lagrangian L = ^{r'^ + r^d"^) -\- ^(r, d) 
with the potential ^(r, ■(?), we linearize the equations with 
respect to x and y. The linearization is valid under the 
assumption that (a) the radial extent of the patch is small 
compared to the scales over which the basic state of the 
galaxy varies significantly and (b) the peculiar motions are 
much smaller than the circular velocities. Note that no as- 
sumption about the circumferential extent of the patch is 
necessary. 



Fig. 1. Amplitude of the potential perturbation, given in ar- 
bitrary units. The initial amplitude of the potential perturba- 
tion at kx = — l.Sfccrit and the minimum at about k^ ~ 1.6fccrit 
correspond to density perturbations of roughly O.OO6S0 and 
O.lSo amplitude, respectively. Since the radial wave number 
increases linearly with time, this also shows the temporal de- 
velopment of the wave. The interval of wave numbers shown 
here, kx — — 1.5 ■ ■ ■ Sfecrit, corresponds to a time interval of 
7 X 10* years. 



This leads to the epicyclic equations of motion 

X = 2noy + AnoAax + f^ (la) 

y^fy- 2xn„, (lb) 

where fx, fy are the acceleration components due to the 
spiral field. Aq denotes Oort's first constant. The pertur- 
bation forces are calculated according to the work of Fuchs 
(1991) where the Boltzmann and Poisson equations were 
used to obtain the perturbation potential in the form 



$(x,y;i) = $o(t)exp(ifc-r) 



(2) 



the following. In the case of free swing amplification one 
has fcj.(i) — k™ + 2kyAf)t for a sinusoidal excitation with 
an initial radial wave number k™ and assuming positive 
ky. While one has to expect that excitations will not be 
sinusoidal in reality, the restriction to a single Fourier com- 
ponent should not be critical due to the linearity assumed 
in the derivation of Eq. (2) and because components with 
wave numbers significantly different from an optimal /c™ 
will either be not amplified strongly or damped out before 
reaching the zone of amplification. 

The integral equation determining the amplitude func- 
tion <I>o(i) derived by Fuchs (1991) is solved numerically. 
Its solution is shown in Fig. 1 for the standard parameters 
of our model. 



fccrit = 27r/(10kpc) 
Q2 = 2 



Ao = 0.015 kms^ipc^i 
Bo = -O.Olkms^^pc"^ 

Ai'VSo = l/(5G7r), 



where k = {kx{t),ky) is the wave vector of the pertur- 
bation. Since swing amplification is most effective in the 
vicinity of ky = 0.5 fccrit, we use this choice throughout 



where fccrit — K^/27rGI]o is Toomre's critical wave number, 
defined with the epicyclic frequency k, the gravitational 
constant G and the unperturbed mass surface density Eo ■ 
Our choice of fccrit gives Eq = 58 Mq pc^^, somewhat more 
than what is found for the Milky Way in the solar neigh- 
bourhood (Kuiken & Gilmore 1991). Q is the Toomrc sta- 
bility parameter, /x'" is the amplitude of the initial den- 
sity perturbation, and Bq denotes Oort's second constant. 
With these parameters one has ilo = 0.025 km s^^pc""'^ 
with a locally slightly falling rotation curve. The only as- 
sumption about the distance to the center of the disk is 
OoTo ^ (7. While the choice of fccrit, ^0, ^0, and Q is 
not meant to accurately model the solar neighborhood, we 
regard it as representative for the outer regions of giant 
spirals. Our results are generic in that they do not qual- 
itatively change under reasonable modifications of these 
parameters. For example, as long as —Bq is comparable 
to Aq, it is unimportant whether the rotation curve is fiat, 
rising or falling. 

Equation (2) describes sinusoidal waves of constant cir- 
cumferential wavelength. Due to the shear of the disk, the 
radial wavenumber kx grows with time. In this way an ini- 
tially leading wave (fc™ < 0) is transformed into a trailing 
one and, as can be seen in Fig. 1, strongly amplified while 
"swinging by" . With increasing kx phase mixing leads to a 
damping of the perturbation as the wavelength of <I> drops 
beyond the typical epicycle size. A single swing amplifi- 
cation event from excitation to decay has a duration of a 
few 10* years. 

In our model the ISM is approximated as a dissipative 
medium of discrete clouds. Their motions are governed by 
equations (1) as long as they do not collide. Collision de- 
tection is performed on a grid of width i?coii, where we 
chose -Rcoii = 50 pc. If there are two or more clouds in 
one cell, we form pairs by a random process and let a 
pair collide if the clouds that make it up are approach- 
ing each other. The collisions treat clouds as "sticky" 
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spheres, i.e. their relative velocities are multiplied by a fac- 
tor (1 — /) after performing an clastic two-body collision 
(Brahic 1977, Schwarz 1981). In such a system the cooling 
rate depends on the collision rate and the inelasticy coef- 
ficient /. To have a clear signature of the coUisonal nature 
of the cloudy medium, we have selected a high collision 
rate (about 0.2stcp"^) which in turn requires a compar- 
atively low / of 0.2 to prevent the system from cooling 
down too quickly. For the implication of this choice, see 
the comprehensive investigation by Jungwiert & Palous 
(1996). However, as long as the total cooling rate is kept 
constant, our model is much less sensitive to the choice 
of / than theirs. With these parameters, an unperturbed 
system of clouds with Gaussian velocity distribution cools 
down exponentially with a characteristic time of about 
200 Myr. Collisions are checked for once every step, where 
one time step equals 3pcskm~^. 

The sticky particles are initially homogeneously dis- 
tributed with a density of 0.0003 pc~^ and have a Gaus- 
sian velocity distribution. We chose a radial velocity dis- 
persion of Uu ~ 4.2kms~^. The velocity dispersion in 
the 2/-direction was set according to the epicyclic ratio 
[2^0 /k)^ resulting in a total initial velocity dispersion of 
(T — 5kms^^. As long as a{t = 0) is low enough (say, 
< lOkms^^), the behaviour of the model is not strongly 
influenced by this choice. 

Velocity dispersions of this order are typical of giant 
molecular clouds (Clemens 1985). Our sticky particles, 
however, are not meant to mimic GMCs themselves. Their 
number is much larger than that of GMCs, they are sta- 
ble, and of course their dynamical behavior can at best be 
regarded as idealized from real clouds. Technically, they 
just implement a dissipativc medium. From the structure 
of the equations of motion, we do not need to make any as- 
sumptions about the masses of the clouds. Also, it is very 
hard to assign radii to the model clouds from the sticky 
particles formalism. One might use i?coii/2 as a gross mea- 
sure of the collision cross section. The details of the model- 
ing do not seem to be critical for our results, since the dis- 
sipativity of the medium does not influence its dynamical 
evolution very significantly except through cooling (see be- 
low). In this way, our results would be largely unchanged 
if, for example, one increased / and lowered the initial 
density accordingly to keep the collision rate constant. 

Introducing two-particle interactions like collisions re- 
quires some care in strictly two-dimensional calculations 
like the present one. Rybicki (1972) showed that in con- 
trast to three-dimensional systems in two-dimensional 
self-gravitating sheets the relaxation time is always of the 
order of the crossing time. Interactions "softer" than grav- 
ity are, however, less affected by this limitation. In fact, 
one can choose Trdax/T'crossing of the ISM freely in our 
model by changing -RcoU and the initial surface density 
since all interactions with impact parameters smaller than 
^coii are equivalent. Checks with 3-D-models showed no 
significant deviations from the evolution described below. 



Fig. 2. The distribution of gas clouds (left) and "backbone" 
stars (right) in a patch of dimensions 20 kpc x 10 kpc, giving 
60000 clouds in the patch. The radial direction is upwards and 
the galactic rotation towards the right. Above each panel the 
step number, the time since the excitation in units of galactic 
rotations (given for illustrative purposes only and computed 
for an assumed ro = 8 kpc) and the kx of the perturbation in 
units of fccrit are indicated. In the frames on the left each dot 
represents a gas cloud and the straight lines mark phases zero 
and ±27r of the shearing wave. The frames on the right side 
are scatter plots of the surface density S(a;, y) obtained from 
inserting (2) into the Poisson equation. To make the structures 
in the stellar disk visible, we enhanced the contrast in the stel- 
lar density by a factor of four. 

3. Results 

3.1. Evolution of a Single Density Perturbation 

In Fig. 2 we show the spatial distributions of clouds and 
stars in a patch of a galactic disk constructed in the way 
described above. We chose a rather large radial extent of 
the patch, bx = 10 kpc, in order to give a clear impression 
of what is going on. The physical relevance of the upper 
and lower edges of the displayed patch may therefore be 
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Fig. 3. Profiles of tlie density p, tfie 
streaming velocities perpendicular and 
parallel to the the wave crest after sub- 
tracting the basic shear, v±_ and iiy, and 
the velocity dispersion, cr, in kms^^ 
across a spiral arm, and the energy 
dissipated in collisions per unit time 
and area, again in arbitrary units. The 
profiles on the left have been taken 
at step 73 and illustrate the state of 
the patch during a formation phase of 
the spiral arms. On the right we illus- 
trate the onset of the dissolution phase 
with profiles taken at step 90. These 
profiles were obtained along lines of 
6kpc length perpendicular to the den- 
sity wave crest. The wave crest is at 
r — 3kpc. 



doubtful, whereas in the region around x = 0, located 
in the center of each frame, the linearization performed 
above is well justified. 

The first frame in the series shows the state of the 
patch about three quarters of a galactic rotation after 
the perturbation has been excited. The perturbation is 
already trailing and well developed in the galaxy's stellar 
"backbone" on the right. The gas has responded to the 
potential perturbation and the amplitude of the density 
perturbation in the gas is comparable to the one in the 
stars. 

This changes about 30Myr later, when the arm- 
interarm density contrast in the stars has just passed its 
maximum of about 1.4:1. As can be seen from Fig. 3 the 
density contrast in the gas is almost 2:1 at this stage. 
The flow pattern is characterized by an inflow towards the 
wave crest and a tangential flow that transports clouds in- 
wards on the outer edge of the arm and outwards on the 
inner edge. Judging from the sinusoidal profiles, we still 
are in the linear regime. These profiles resemble analo- 
gous results from simulations with discrete clouds in quasi- 
stationary density waves (the results of Roberts (1992) are 
particularly well suited for comparison). The main differ- 
ence is that shearing density waves are essentially at rest 
with respect to the disk material around them and thus 
produce structures symmetric to the wave crest, whereas 
stationary density waves outside the corotation zone have 
an upstream and downstream side and thus lead to asym- 
metric profiles. 

While the perturbation amplitude in the stars declines 
from now on (cf. Fig. 1), the spiral arm in the gas continues 
to grow for another 50 Myr, as seen in the third frame in 
Fig. 2. With this maximal density contrast in the gas. 



the formation phase of the spiral arm comes to an end. 
The following dissolution phase is marked by a gradual 
turnover from infiow to outflow with respect to the wave 
crest. In the velocity component perpendicular to the wave 
crest, v±, a double- wave pattern evolves (cf. Fig. 3). 

As can be seen in the fourth frame in Fig. 2, the spiral 
arms widen during the dissolution phase, and their edges 
sharpen. This shift of activity towards the edges is dis- 
cernible in the cr-proflle shown in Fig. 3, right side, where 
one already has maxima at the arm's edges. At this stage 
hardly any perturbation is visible in the stellar disk. 

In the last frame shown in Fig. 2, taken about half 
an epicyclic period after the strongest expression of the 
pattern, the gas clouds once more gather into a density 
concentration, this time at phases ±(2j -I- l)7r, j S IN, of 
the original density perturbation. Apart from the shift in 
phase and the smaller wavelength, proflles in p, vj_, and 
W|| across this "echo" arm closely resemble those given in 
Fig. 3, left column. However, it has no visible counterpart 
in the stellar disk. Another half epicyclic period later one 
again finds an echo, this time with maxima at the positions 
of those of the original density perturbation. These echoes 
are only slowly damped out. 



3.2. Physical Interpretation 

To get an understanding of the dynamics of the patch, it 
is advantageous to introduce the crest distance in units of 
the perturbation wave length, r(t) = rfe/(27r). In f unper- 
turbed circular orbits appear as parallels to the x-axis, and 
the acceleration due to the perturbation is always parallel 
to the measuring rod. The drawback is that the length 
unit varies with time, so that orbits of constant nonzero 
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Fig. 5. Time evolution of crest distances. Here 
we show the distances of 15 collisionless parti- 
cles starting with zero peculiar velocity (solid 
lines) in units of the time-dependent wave- 
length of the potential perturbation. The long- 
dashed line connects the epochs of the turnover 
of the particles' motion, and the vertical line 
marks the time at which the profiles on the 
right hand side of Fig. 3 were taken. The short- 
dashed lines indicate crest distances of two gas 
clouds from our reference simulation for com- 
parison. Note that A(f) varies between 20kpc 
and 2 kpc in the time interval shown in this di- 
agram, hence the amplitude of the oscillation 
is no measure for the sizes of the epicycles. 



epicyclic amplitude appear to have growing amplitude in 
f when kx > 0. 

The evolution of r{t) for 15 particles moving in our 
reference model without collisions is shown in Fig. 5. The 
particles start on circular orbits and are accelerated to- 
wards the wave crest for the first 200pcskm~^, indi- 
cated by a negative curvature in f(t). At about the time 
when the stellar arm starts to dissolve {k^ « l.T/ccrit by 
Fig. 1, corresponding to i « 210pcskm^^), the curvature 
changes sign. The net force on a particle must therefore 
be directed away from the wave crest, although $o(0 is 
still negative. This is because the Coriolis force present 
in the equations of motion (1) is now stronger than the 
perturbation force. While the latter decays, the Coriolis 
force finally turns around the motion. Since the relative 
strengths of perturbation and Coriolis forces vary with 
the crest distance, the epoch of this turnover also depends 
on f. This delay in turnover time causes the double wave 
shown in Fig. 3 on the right side. 

While it is possible to gain a qualitative understanding 
of the variation of turnover times in terms of the ratio of 
perturbation versus Coriolis force, the problem is better 
understood in terms of a perturbation of the epicyclic fre- 
quency, i.e. by letting the angular frequency of a supposed 
periodic motion in a shearing reference frame vary with 
the distance to the wave crest. For a simplified version of 
the present potential, we derive 



K?{r) 



ka 
1 coslkr] 



(3) 



in the appendix for the perturbed epicyclic frequency k 
where a is the amplitude of the perturbation — which is 
negative during the formation phase — , k is its wave num- 
ber and r denotes the distance from the wave crest. Thus, 
the double wave is a consequence of k increasing towards 
the wave crest. As can be seen from Fig. 4, the rough es- 
timate Eq. (3) already is quite a good approximation to 
the more complex situation in our simulation. 



Fig. 4. Perturbed epicyclic frequencies k for two distances 
from the wave crest as a function of time. The solid lines 
are computed from Eq. (3). The dashed lines are averages for 
the epicyclic frequencies over the numerically integrated or- 
bits of 250 particles distributed evenly on an epicycle with u = 
10 km s~^ . For the upper curve, the guiding center of the epicy- 
cle lies at a:: = 200 kpc, for the lower curve at a; = 1500 kpc. 



Fig. 5 also shows how the dissolution phase turns over 
into the formation of the echo arm as the clouds pro- 
ceed on their epicycles. After t — 300 Myr {k^ ~ 3A:crit5 
cf. Fig. 1) the potential perturbation is practically zero, 
and therefore the echoes simply are kinematical spiral 
arms (Kalnajs 1971). They do not appear in the stellar 
disk because of the much higher velocity dispersion of the 
stars. For our choice of ky = 0.5 fccrit, the size of the stellar 
epicycles gets comparable to the perturbation wavelength 
shortly after the maximal amplification, and phase mix- 
ing leads to an exponential damping. Thus the only sign 
of kinematical spiral arms in the backbone stars is a small 
overshoot of $0 around kx = 4A:crit (cf. Fig. 1). 

In contrast, the echo arm in the gas has about the am- 
plitude of the primary spiral arm. This is mainly because 
of the smaller velocity dispersion of the gas. When phase 
mixing becomes important in the stars, the wave length 
of the perturbation is still much larger than the typical 
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size of an epicycle of a cloud. Only at a wavelength of 
about 500 pc (t > 900 Myr) does phase mixing dominate 
the damping of the perturbation in the gas. Until then, 
the perturbation is mainly damped by collisions that are 
not very efficient in damping away larger-scale structures. 
Thus, even when radial variations of k weaken the kine- 
matical spiral arm, at least the first echo can be expected 
to be reasonably strong. One might speculate that echoes 
of this kind have some relation to "gaseous interarms", 
spiral arms in the ISM without corresponding activity in 
the stellar disk (e.g.. Block et al (1994), although their 
observations only cover spirals with much higher density 
amplitudes than we consider here). 

The discussion of the dynamical evolution of the patch 
given so far did not make any reference to the dissipativity 
of the cloudy medium. In fact, turning collisions off does 
not strongly change the evolution of the density pertur- 
bation in the ISM. The impact of the dissipativity on the 
dynamics is small because the streaming motions induced 
by the perturbation are quite smooth, so that the dissi- 
pativity of the medium mainly appears as viscosity that 
can be disregarded on short time scales. On longer time 
scales, however, the dissipativity of the cloudy medium 
keeps its velocity dispersion low, whereas a single swing 
amplification event roughly doubles the temperature of a 
stellar system with cr(t = 0) = 5kms^^. 

Spiral arms are traced by young stars and HII- 
regions. Although the question has been under some de- 
bate (Elmegreen & Elmegreen 1986), it appears that den- 
sity waves do modulate either the initial mass function 
or the star formation rate nonlinearly (Cepa & Beckman 
1990). Kwan & Valdes (1983) have argued that massive 
star formation is triggered by collisions between diffuse 
molecular clouds. Following this line one might expect the 
energy dissipated in collisions per unit time and area, 77, to 
give a coarse measure for the efficiency of massive star for- 
mation. Though this is a naive view of how star formation 
proceeds in molecular clouds, the compression induced by 
the development of spiral arms certainly plays a role in 
triggering the collapse of molecular clouds and hence star 
formation. As can be seen from Fig. 3 we find contrasts in 
?7 of up to 1 : 10 from arm to inter- arm regions since j] is 
larger in the arms. 

3.3. Consecutive Swing Amplification Events 

Certainly the processes investigated here will not take 
place only once for each galaxy; instead, we expect ex- 
citations to recur in a more or less stochastic manner. To 
examine the behavior of the patch under repeated per- 
turbations, we need to make some assumptions regard- 
ing the frequency and the amplitude spectrum of the am- 
plification seeds. In our model, we are rather limited in 
the first of these choices because the crests of the den- 
sity waves must always connect the centers of adjacent 
patches to keep the perturbation potential continuous 



across patch borders. Consequently, the swing amplifi- 
cation events must follow one another with a period of 
by/{2Aobx) (or a multiple thereof), where the circumfer- 
ential extent of the patch is already fixed to by — 2TT/ky by 
the continuity condition across circumferentially adjacent 
patches. 

In contrast, we expect swing amplification events to 
follow one another in a more or less random fashion. This 
should be the case for seeds provided by Poisson noise or 
massive perturbers of the disk (Toomre 1990). Even recur- 
rent swing amplification due to mode coupling as stud- 
ied by Fuchs (1991) will have a stochastic nature. Thus 
strictly periodic excitations are probably grossly unphys- 
ical. However, as neither the global kinematics nor the 
overall morphology of the patch strongly depends on the 
excitation period as long as one avoids resonances with 
the epicyclic period, one can still gain some insight in the 
long term evolution of the patch. For simplicity, we use 
the same amplitude for all seeds. 

The general appearance of the patch subjected to 
such repeated perturbations qualitatively resembles the 
sequence shown in Fig. 2 for a single swing amplifica- 
tion event. However, the spiral arms are more ragged with 
crossings and bridges, because the potential perturbations 
act on a disk that is already structured both in density and 
velocity. For excitation periods used here, only one swing 
amplification event causes noticable perturbation forces at 
any moment. This corresponds to the fact that the lifetime 
of any perturbation will be much shorter in the backbone 
stars with their high velocity dispersion than in the cool 
cloudy medium. One consequence of this is that the spiral 
arm in the backbone stars will be much smoother than in 
the cloudy medium with no crossings and bridges appear- 
ing. Again, this invites comparison with the observations 
of Block et al (1994), who see much smoother arm shapes 
in the infrared than in the optical. 

In longer-term simulations the gas reaches a quasi- 
equilibrium state after the first few excitations. Then, the 
global velocity dispersion including streaming motions os- 
cillates by about Ikms^^ around a mean value of about 
a w 6kms^^ for our parameters. This mean value in- 
creases slightly with decreasing Tg over a range of exci- 
tation periods from 150pcskm~^ to 400 pes km" ^ (ex- 
cluding resonances with k). This equilibrium is remark- 
ably independent of the collision parameters. Going from 
a mean collision rate of 0.14step~^ to 0.3step"^ one finds 
the temporal mean of the velocity dispersion decreases 
by less than 0.25 km s"^, where most of the difference is 
due to the systems with high collision rates cooling down 
more quickly from a peak velocity dispersion of about 
6.5 km s"^ that is almost independent of the collision rate. 
The reason for this behaviour is that in the later evolu- 
tion a is dominated by the thcrmalization of streaming 
motions. These are to zeroeth order independent of the 
velocity dispersion before the swing amplification event. 
Thus, cool systems are heated to the about the same tem- 
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perature as warmer ones. As long as the system is able 
to cool down between two swing amplification events, and 
provided that the streaming motions are not destroyed 
by too frequent collisions, the sole difference between the 
cooler and warmer systems lies in the the amplitude of 
the temperature oscillations. Only when collisions become 
unimportant — at collision rates of, say, 0.05 step^^ — , does 
this reasoning break down and the mean velocity disper- 
sion rises to values of the order of lOkms^^. The velocity 
dispersion of a collection of coUisionless test particles in 
the perturbed disk rises to about 20kms^^ within 3 Gyr, 
which is still much less than the 45 kms^^ of a stellar disk 
with the values of Q and Sq we assume. 



Fig. 6. Profiles of p and v±_ analogous to the ones shown in 
Fig. 3, this time taken after 1.5 Gyr of evolution under consec- 
utive swing amplification events. 



The profiles in Fig. 6 show that by constructive in- 
terference between two swing amplification events a pre- 
structured ISM can show a much more vigorous response 
to a perturbation than the homogeneous patch exam- 
ined above. Of course, destructive interference also occurs. 
Finding profiles like the one shown in Fig. 6 should be 
quite within reach of current observational techniques. 

One also finds that the double- wave pattern can persist 
for up to 1/3 of the lifetime of a spiral arm, depending on 
how the excitations follow one another. 



3.4. Self gravity 

Jog (1992) predicts from her two-fluids model that the 
oscillation frequency for spiral arms in the gas is much 
lower than for the stars from a self-consistent two-fluid 
model. This is not the case in our cloudy medium. In fact, 
kinematical spiral arms in an ISM consisting of massless 
clouds will oscillate with the epicyclic frequency. The ques- 
tion arising naturally is whether the oscillation frequencies 
in gas and stars are so similar because of the neglect of 
self gravity in the ISM. 

The investigation of the effects of self gravity could 
be implemented by a straightforward iV-body code for 



our discrete clouds. However, considering Fig. 2 one no- 
tices that the density evolution in the gas has the form 
of a shearing wave that closely follows the evolution of 
the stellar perturbation in wavelength and inclination. 
Thus the additional perturbation potential arising from 
the gas clouds can be taken into account by assuming 
that the surface density of the gas sheet is modulated like 
/^gas oc exp(ifc ■ r) with k = (k™ + 2AQkyt, ky). 

Assuming a linear response the Poisson equation yields 
the additional potential perturbation arising from self 
gravity 



$scif (a;, y; t) = K 



|fc'"| 6S0 

k{t)\ yU'" 



$o(fci")cos(fc. 



(4) 



where fc'" denotes the wave vector of the initial pertur- 
bation, e — /igas/Mstars and K is the Fourier amplitude of 
the perturbation in the gas on cos(fcr). In practice, K is 
determined by a least-squares fit of K cos(fc • r) + /I to a 
density cut perpendicular to the wave crest. 

The assumption of a linear response is certainly valid 
in the early stages of the evolution, while its justification 
is somewhat doubtful for the later stages of the dissolution 
phase when the edges of the arm get very sharp. Still, one 
may expect results that allow a qualitative judgment on 
the importance of self gravity. 



Fig. 7. Time evolution of the density contrast of the spiral arm 
for three values of the gas-to-stars mass ratio. For a low gas 
mass fraction, the effect of self gravity is negligible. 

Figure 7 shows that neglecting self gravity is well jus- 
tified in our model for gas mass fractions similar to those 
found in the solar neighborhood; the oscillation frequency 
of the arm in the gas does not change much for e = 0.1. 
For late-type spirals with a high gas content, however, self 
gravity alters the behavior significantly and our results 
should be applied with care. In particular, looking at pro- 
files analogous to Fig. 3 one finds that a higher gas mass 
fraction will lead to softer edges and a less pronounced 
double wave during the dissolution phase. 

Of course, this scheme, being designed to assess the 
consequences of neglecting self gravity on the large scale 
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evolution in our two dimensional sheet, cannot describe 
many effects certainly important to disk evolution. In par- 
ticular, Toomre (1990) pointed out that an A^-body sim- 
ulation in a patch similar to ours rapidly develops GMC- 
like complexes that might act as seeds for swing amplifica- 
tion. One can expect these to occur as well in the cloudy 
medium investigated here if it were fully self gravitating. 
However, an examination of these processes, possibly lead- 
ing to self-triggered swing amplification with gas coupling 
back on the stars, is beyond the scope of this paper. 



4. Conclusion 

In this work we used a model consisting of a stellar con- 
tinuum and gas clouds behaving like sticky particles to 
examine the behavior of gas in transient density waves. In 
contrast to the prior work of Jog (1992) that employed a 
gas-dynamical approach we find that the evolution of the 
density perturbation in the gas remains closely coupled to 
the one in the stars even when self gravity of the ISM is 
taken into account. 

In a first phase, the ISM responds to the growing po- 
tential perturbation by flowing towards the wave crest. 
When the potential perturbation already decreases, the 
amplitude in the gas continues to grow until the dissolu- 
tion of the gaseous arm begins due to Coriolis forces. In 
this process, a variation of a cloud's epicyclic frequency 
with its distance from the wave crest leads to a double 
wave in the velocity perpendicular to the arm, with par- 
ticles near the wave crest already flowing outward while 
those further out are still moving in. We regard this as a 
strong signature for transient density waves that might be 
within reach of current observational techniques. Finally, 
the gaseous arms continue to evolve as kinematical spiral 
arms until they are finally damped out by phase mixing. 

For the longer-term evolution we propose that swing 
amplification events follow one another in a more or less 
random fashion. We find that for a wide range of excitation 
frequencies spiral arms continue to form and dissolve. By 
interference between density waves a ragged morphology 
results. The dissipativity of the gas — unimportant for the 
short-term dynamics — keeps the velocity dispersion of the 
clouds low and thus their responsiveness high. 

Given the stochastic nature of the processes described 
here, one cannot predict the appearance of a galaxy based 
on a knowledge of its basic state using the methods ap- 
plied in this work. However, since our results are quite 
insensitive to the choice of parameters of the underlying 
model, we regard them as fairly generic for transient den- 
sity waves. Thus, finding a double wave in radial velocity 
in the line of sight across a spiral arm for some component 
of low velocity dispersion (e.g., gas or OB-associations) 
would make a strong point for the existence of transient 
spiral arms. 
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Appendix 

In the special case of a weak and purely radial per- 
turbation, {fx,fy) = (asin(A:a;),0), it is easy to show the 
dependence of the local oscillation frequency on the posi- 
tion of the particle within the density wave. While this is 
quite different from the potential arising from a shearing 
density wave, one can readily verify by numerical integra- 
tion of the orbits that the oscillation frequency of mass 
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points in a potential of the type of Eq. (2) roughly fol- 
lows the behaviour predicted here with minor deviations 
resulting from the time-dependence of the wavelength of 
the perturbation (Fig. 4). 

The assumption of a purely radial perturbation allows 
a closed integration of Eq. (lb). Inserting the result into 
Eq. (la) yields 



where Jq and Ji denote the Besscl functions of order 
and 1, respectively. By comparing the 'leading' terms pro- 
portional to sn\{kt) one obtains 



«2 



1 



2a 



Ji{kxi) cos(fca:;o~ 



(A3) 



X = 2ilvo — K X + a sm{kx) 



(Al) 



A straightforward application of (A3) to the pertur- 
bation used in our simulations — ^justified a posteriori by 
Fig. 4 — can be done by identifying xq with the distance of 
the guiding center of a particle from the wave crest. Since 
the epicycle sizes are small compared to the wavelegth 
27r/fc of the perturbation, one may furthermore approx- 
imate Ji{kxi)/xi « k/2. This leads to Eq. (3). These 
approximations obviously break down for epicycle sizes of 
the order of a typical stellar epicycle size where the am- 
plitude of the oscillation in k has to be zero. 



where vo is an integration constant. To the author's knowl- 
edge, a closed solution of this differential equation does 
not exist. We therefore chose to investigate its solutions 
using the method of harmonic balance (e.g., Bogoliubov & 
Mitropolski 1961). We assume a quasi-harmonic solution 
of the form 

x{t) — xa + xi sin(Ki), 

where R = k{xo,xi). Inserting this into Eq. (Al) and 
Fourier analyzing the nonlinear term leads to 

—K Xi sm{kt) — K xo — 2floVQ + asm(kxx)Jo{kxi) 

r, This article was processed by the author using Springer- Verlag 

+ (k^ cos(fca;o)Ji(fca:;i))a;i sin(Kt)(A2) TJjX A&A macro package 1992. 
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